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We discuss the discrete spectrum of the Hamiltonian describing a two- 
dimensional quantum particle interacting with an infinite family of 
point interactions. We suppose that the latter are arranged into a star- 
shaped graph with N arms and a fixed spacing between the interaction 
sites. We prove that the essential spectrum of this system is the same 
as that of the infinite straight "polymer" , but in addition there are 
isolated eigenvalues unless N = 2 and the graph is a straight line. We 
also show that the system has many strongly bound states if at least 
one of the angles between the star arms is small enough. Examples of 
eigenfunctions and eigenvalues are computed numerically. 



1 Introduction 

Graph-type systems are used in quantum mechanics for a long time 
but only in the last decade they became a subject of an intense interest - 
cf. |[KS| | and references therein. Among various graph geometries, star graphs 
were investigated from different point of view. Recall, for instance, a natural 
generalization of the weak- coupling analysis for one- dimensional Schrodinger 
operators ||E1||, signatures of quantum chaos found recently in stars with finite 



nonequal arms [BBK], etc. 
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From the mathematical point of view Schrodinger operators on graphs are 
easy to deal with, because they represent systems of Sturm-Liouville ODE's 
coupled through boundary conditions at the graph vertices. This is due to 
the assumption that the configuration space of the system is just the graph. 
From the physical point of view, this is certainly an idealization. One of 
the most common applications of graph models is a description of various 
mesoscopic systems like quantum wires, arrays of quantum dots, etc. In 
reality their boundaries are finite potential steps, and therefore the particle 
can move away from the prescribed area, even if not too far because the 
exterior of such a graph is a classically forbidden region. 

There are various ways how to model such "leaky" graphs. One can 
use Schrodinger operator with a Dirac measure potential supported by the 
graph - see [|BT| , [Ell and references therein. Here we consider another, in 
a sense more singular model where the graph is represented by a family of 
two-dimensional point interactions. Its advantage is that such a model is solv- 
able because (the discrete part of) the spectral analysis is reduced essentially 
to an algebraic problem. Two-dimensional point-interaction Hamiltonians 
were studied by various authors - references can be found in the monograph 




Nevertheless, relations between spectral properties of such oper- 



ators and the geometry of the set of point-interaction sites did not attract 
much attention. Here we are going to fill this gap partly by discussing an 
example of a point-interaction star graph. 

The model is described in the following section. Next, in Section 3, we 
show that the essential spectrum is given by the structure of each graph arm 
at large distances and thus it coincides with that of an infinite straight "poly- 
mer" |[AGHH| , Sec. III. 4]. More surprising is the fact that a star graph has a 
nonempty discrete spectrum, with the exception of the trivial case when the 
graph is a straight line. This is proved in Section 4 where we also show that 
there are geometries which give rise to numerous strongly bound states. In 
the final section we present numerically computed examples showing eigen- 
values and eigenfunctions for various graph configurations. Of course, the 
discrete spectrum is not the only interesting aspect of these Hamiltonians. 
One can ask about the scattering, perturbations coming either from changes 
in the geometry or from external fields, etc. We leave this questions to a 
future publication. 
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2 Formulation of the problem 



For a given integer N > 2, consider an (iV — l)-tuple of positive numbers 
(3 := . . . , Pn-i) such that J^fji (3j < 2n and denote i?j := ^f =1 A an d 
$o := 0. Then one can define the set 

JV-l 

Y= \J{{nl cos(^), nl sin(^))} neN U (0, 0), 

where I > is a given distance which has the meaning of the spacing of 
points at each "arm" of Y. 

The object of our study is a two-dimensional Hamiltonian, which we 
denote as Hjsr(a, (3), with a family of point interactions supported by the set 
Y having the same "coupling constant" a. The point interactions are at that 
defined in the standard way ||AGHH|| by means of the generalized boundary 
values, 

1 

L (i/j,a) := lim — — — -rip(x), 
[x-3|-»o In \x — a\ 

Li(if),a) := lim [if)(x) — L (i[), a) In \x — a\] . 

|af — a| — >0 

Due to its point character, the Hamiltonian acts as free away of the inter- 
action support, (H N (a, f3)ip)(x) = (— A-0)(x) for x $ Y, and its domain 
consists of all functions ip G W 2,2 (M. 2 \ Y) which satisfy the conditions 

la C0, a) + 27tq;L o ('0, a) = 

at any point a from the set Y. Since the particle mass plays no role in the 
following, we choose the units in such a way that 2m = 1. 



3 The essential spectrum 

Consider first the essential spectrum of Hn{oc, (3). It is well known for the 
so-called straight polymer, i.e. H2(a,ir), which is discussed in [ |AGHH| , 
Sec. III. 4]. In this particular case the spectrum is purely absolutely con- 
tinuous and has at most one gap. Specifically, it equals [E ,Ei] U [0, oo), 
where E < E\ < 0, for the coupling stronger than a critical value, a < ay, 
while in the opposite case the two bands overlap, E\ > 0, and the spectrum 
covers the interval [E , oo). The values E ,Ei, and ay are given as implicit 
functions of the parameters a and I. 
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Proposition 3.1 The relation inf a ess (H N (a, f3)) = inf a(H 2 (a, tt)) holds 
for any (3 and N . 

Proof: The easy part is to check the inclusion a ess (H N (a, (3)) D a(H 2 (a, it)). 
Given an arbitrary A G a(H 2 (a,ir)) we construct a sequence {ip n }'^ = i with 
i>n(x) = jn(x)4>x(x), where (fix is a generalized eigenfunction of H 2 (a, n) with 
the energy A and j n G Cq°(M. 2 ) are mollifier functions to be specified. If 
supp ] n intersects just one arm of Y, we have 



Since the functions (fix, V(fix are bounded, it is sufficient to take j G C^°(IR 2 ) 
with 1 1 j 1 1 = 1 and define 



for a suitable sequence {x n } C M 2 ; the latter can be always chosen in such a 
way that each j n intersect with a single arm of Y. Using || Vj ra || = n _1 / 2 || Vj|| 
and ||Aj n || = n _1 ||Aj||, we conclude that (H N (a,(3) — \)ip n — > strongly as 
n — > oo, i.e. that A G a(H N (a, (3)). We can even choose {x n } so that ip n 
have disjoint supports forming thus a Weyl sequence, but it is not needed, 
because a(H 2 (a, 7r)) consists of one or two intervals and A belongs therefore 
to the essential spectrum of Hn{ol,P). 

To prove the inequality inf a ess (HN(o:, f3)) > inf a(H 2 (a, 7r)) we employ 
the Neumann bracketing. We decompose the plane into a union 



where Sj is a half-strip centered at the line {x G M 2 : argx = -dj} of the 
width d, Wj is a wedge of angle (3j + \ between two half-strips Sj and Sj+i, and 
finally, P is the remaining polygon containing the center part of the "star". 
Introducing the Neumann boundary conditions at the boundaries, we obtain 
a lower bound to H^(a,P). This new operator H is equal to a direct sum 
of Neumann Laplacian corresponding to the said decomposition. Since each 
wedge part have a purely continuous spectrum equal to R + and the 
polygon part H P has a purely discrete spectrum, the half-strip parts H Sj are 
crucial for the threshold of the essential spectrum of H. 



(H N (a, (3) - \)ij n = (H 2 (a, tt) - \)^ n = -2V0 A • Vj n - A Aj,, 





(3.1) 
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We can choose the boundaries such that the distance between the trans- 
verse boundary of the half-strip and the first point interaction is equal to 
1/2. The spectrum of Hs 3 on this half-strip is the same as the symmetric 
part of spectrum of a Neumann Laplacian Hd on a "two-sided" strip of width 
d, hence the threshold of cr e < jS (ifs j ) coincides with that of a ess (Hd). 

Following the standard Floquet-Bloch procedure - see | AGHH , Sec. III. 3] 
- we can pass from Hd to an unitarily equivalent operator which decomposes 
into a direct integral, 



UH d U~ l = ^r H d {6)d6 
171 Jee\-TT/i,TT/i) 



Hr/WQ 

where Hd(6) is a point-interaction Hamiltonian in L 2 ([0,Z] x [— d/2,d/2]) 
which satisfies the Bloch boundary conditions, 



4<{l-,y) = e m ii0+,y) 



OX ox 



for y e [—d/2, d/2}. The position of the point interaction is chosen as (a, 0). 
Then it is easy to write the corresponding free resolvent kernel with one 
variable fixed at that point, 



oo oo 



Gq(x, a; 9, z) 



■I 2irm 



+e){x-a) 



m=—oo n=0 



X COS 7171 



y + d/2 
d 



cos titc 



_0 + d/2 
1 



Using the formula [BMP, 5.4.5.1] one can evaluate the inner series getting 



1 oo 

G d (x,a;9,z) = \ ^ e <^^)(— ) 



m=— oo 
1 



ra\' J i 



cosh((rf - |y|)x tn (6>, z)) + cosh(yx m (9, z)) 
2 smh(dx m (8, z)) 



where x m (0,z) = ^/(^p + ^) 2 ^. 

To compute the generalized boundary values L and L 1; and from them 
the eigenvalues of Hd(0), we follow the procedure from |[KGST|| . The coeffi- 
cient at the singularity does not depend on the shape of the region [Ti|, i.e. we 
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have Lo(^, a) = —l/2m/j(a). The value L\ is expressed by means of the regu- 
larized Green's function, £(e; 9, z) := rim|£_g|_> (Go(x, a; 9, z) +ln \x — a|/27r), 
where we introduced the e = l/d with a later purpose in mind; to compute 
it we replace the term In \x — a\ by its Taylor series and perform the limit 
x — > a under the series. 

Recall that we are interested in the lowest eigenvalue of Hd(0), and that 



due to general principles ||Wej , Sec, 8.3] a single point interaction gives rise 
to at most one eigenvalue in each gap for a fixed 9, and this is given as a 
solution to the implicit equation a — £(e, 9, z). Putting 9 = 0, we have for z 
in the lowest gap 



— oo 
oo 



+ 



l 



cosh(e" 



m=l 



2W 
cosh(e _1 



—z sinh(£ z) 
x m (0,z)) + l Z 



x m (0, z) sinh(e- 1 x m (0, z)) 



2nm 



So far we made no assumption about d. It is obvious that the decom- 
position ( p.ip can be chosen in such a way that d is an arbitrarily large 
number. For a large d, i.e. small e, the first two terms of Taylor series for 
the eigenvalue read 



z\e, a 



dz 



e + 0(e 2 



where the first derivative is easily computed by means of the implicit-function 
theorem, 



dz(e, a) 



de 



yc 



(^P) 2 -^(o,«) 



y c 



(2zp) 2 _ z(0 ,«) 



3 ' 
2 



it is well defined and negative for z(0, a) < 0. 

Hence z(e, a) approaches z(0, a) from below as e — > 0; it remains to prove 
that the limit value is the threshold E . To this aim, we apply the Floquet- 
Bloch decomposition to the operator H2(a,n). The formulae for Green's 
function and the ^-function at 9 = change to 



Gnix, a:9,z) 



oo 



21 2— 



-x m (e,z)\y-b\ 



x„ 
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and 

1 / " 

x m (0, z) 2nm\ ' 

respectively. It is obvious that the solution to the equation a — equals 
z(0, a). Summing the argument, we found that inf a(Hd) — > inf a(H 2 (a, ir)) 
from below as d — > oo, and therefore inf a ess (H N (a, (3)) > inf a(H 2 (a, n)) — r\ 
holds for any 77 > 0, which concludes the proof. I 

If a > ay the above proof shows that a ess {H^{Q.^ 0)) = o~(H2(a,n)), 
while in the opposite case one should check also that the two spectra have 
the same gap. Since the coincidence of the two spectra is not important in 
the following, we are not going to discuss this question here. 

4 The discrete spectrum 

Our main claim in this paper is that the geometry of the star graph gives 
rise to a nontrivial discrete spectrum, and that for some configurations there 
are many strongly bound states. We shall state the result as follows: 

Theorem 4.1 (a) adi SC (H N (a, (3)) 7^ unless N = 2 and (3 = n. 

(b) Let N, I, and a be fixed. For any positive integer n and c G M. one 
can choose the graph geometry (making at least one of the angles (3j small 
enough) in such a way that the number of eigenvalues of H N (a, (3) (counting 
multiplicity) below c is not less than n. 

Proof: As we mentioned above, the straight polymer has empty discrete 
spectrum. The existence of at least one eigenvalue below the threshold Eq 
for N = 2 and (3 7^ 7r has been established in |[E2|| . The part (a) then follows 
from a simple auxiliary result. Let o~disc(H) = {Ej : E\ < E 2 ■ ■ ■ < -EW } 
with iV finite or infinite be the discrete part of the spectrum of a self-adjoint 
operator H. We will say that o-& m (H') < a disc (H) for another self-adjoint 
operator H' if ^a disc {H') > #o- disc (H) and E\ < E j for all j = 1, ... , N. We 
claim that adding an arm to the graph pushes the discrete spectrum down, 
adding possibly other eigenvalues on the top of the shifted eigenvalue set. 

Lemma 4.2 crdisc(-ftjv(a, (3)) > (Xdisc(-£^v+i(tt, pj) holds for any N and angle 
sequence (3 = (A,..., fa. x , pf\$f\ (3 j+1} f3 N ^) with ^ + /3f = (3 r 
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The proof would be easy if the two operators allowed a comparison in the 
form sense and the minimax principle could be used. It is not the case, but 
one can get the lemma by induction from the following result. As in [ |AGHH|| 
we denote by H a y the Hamiltonian with point interactions supported by 
an arbitrary set Y = {%■}; we suppose here that all of them have the same 
coupling constant a. 

Lemma 4.3 o- disc (H aX >) < a disc (H aX ) for any Y' = Y U {y'} with y' g Y. 

Proof: By[ |AGhfH| , Thm. 1.5.5] H aY can be approximated in the norm- 
resolvent sense by a family of Schrodinger operators H e a Y = —A + V e with 
squeezed potentials supported in the vicinity of the points of Y, and the same 
is true for H aY i. Each part of the potential can be chosen non-positive, and 
the parts corresponding to the points of Y may be the same for both approx- 
imating operators, in which case we have H* Y , < H e a Y in the form sense, 
or even in the operator one if V £ and V' e are regular enough. By minimax 
principle we infer that a d i SC (H^ y) — °~disc(Ha y) holds for any e > 0, and the 
relation persists in the limit e — > in view of the norm-resolvent convergence. 
I 



Proof of Theorem \4.1[ continued: By Lemma fO] it is sufficient to prove the 
part (b) for N = 2. We choose then j3n < 1/2 which makes it possible to draw 
n circles of radius R < 1/2 centered at the points (jl, 0), j — 1, . . . , n. Each of 
them contains exactly two point interaction, those placed at (jl, 0) = Si and 
(jl cos (3, jl sin pi) = a 2 . Their distance is therefore a = |ai— a 2 \ = 2jl sin(/3/2). 
By imposing Dirichlet boundary condition at the circle perimeters, we obtain 
an operator estimating H^(a,f3) from above. 

The proof is now reduced to the spectral problem of the Hamiltonian 
H(a, a) with two point interaction in a circle with the Dirichlet boundary; 
it is sufficient to show that to a given c there is ao > such that H(a, a) 
has an eigenvalue < c for each a G (0, ao). This operator has at most two 
eigenvalues which are solutions of the following implicit equation, 

det A(a,a 1 ,a 2 ; z) = (4.2) 

with 

A ij (a,a 1 ,a 2 ;z) := 5ij(a - ^(a,; z)) - (1 - 5 ij )G (a i ,a j ; z) , i,j = 1,2, 
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where Go is the integral kernel of the resolvent (H(a, a) — z)^ 1 and £(a.j; z) 
is the regularized Green's function obtained by removing the logarithmic 
singularity at a^. Due to the rotational symmetry of the region we have 

G {ax,x;z) = ^- ( K (x\x - a x \) - 7 H^-ai|) ] , ( 4 - 3 ) 

where x = \/— z. Since it is sufficient to consider c < £7q < we may suppose 
that z is negative. The first ^-function is easy to compute, 

= -^(D-h-.^j. (4.4) 

The second one is more difficult, because to express it one would need to 
replace the formula (|4.3|) by the Green function with a general pair of argu- 
ments. Instead we employ a simple Dirichlet bracketing argument similar to 
that used in flENl. 



Consider three Hamiltonians with a single point interaction placed at a^'- 
with no restriction in the whole plane, Ho, with the Dirichlet condition 
at the circle with center a\ and radius R (the same as for H(a,a)), and 
finally H + , with the additional Dirichlet condition at a circle with center a 2 
and a radius R 1 < R — a. These operators satisfy obviously the inequalities 
H < Hq < H + ; since we are interested in comparing the negative spectra, 
only the interior parts of the last two operators have to be considered. The 
inequalities between the ground-state eigenvalues of the three operators imply 
inequalities for corresponding ^-functions, 

£+(a 2 ; z) < £(a 2 ; z) < £_(a 2 ; z). 



Both £± are known, one from | [AGHH| ], the other from ( [4.4] ) with changed 
parameters. In this way we get 



i{a 2 -z) = J-^(l)-ln|-C(x: 



2tt V r v ' 2 

with the "error term" satisfying < C(x) < jfr^Tj ■ 

One can choose xo > in such a way that 2ira — ip{l) + ln(xo/2) is 
positive. Using the monotonicity of u i— > ^(u) anc ^ ^he aDove inequalities for 
C(x), we get from the condition (^4.2| ) for x > x the estimate 



x 



2m _ m + hl _ + j > ^ o(xo) _ 
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We employ further the behavior of the modified Bessel functions 1$ and K 



as xa — > 0, see [BMP , 9.6.12-13]. For small enough a we can thus choose 
the positive square root of the r.h.s. and the condition ( |4.2| ) has a solution 
satisfying 

X > A exp L (l ) 7J-Q; f^fH (1 + 0(a)). 



Since the corresponding eigenvalue is — x 2 and all the a in the estimating 
operators can be made simultaneously small by choosing (3 small enough, 
the proof is finished. 1 

Notice that for small a the estimating operators have only one eigenvalue. 
Considering two point interaction in the whole plane, we see that the estimate 
is reasonably good: we have 

x < -p e m ~™ (1 + 0(a)). 



5 Numerical results 
5.1 The Method 



By [|AGHH| , Thm. III. 4.1] the Hamiltonian Hn(o!,(3) can be approximated in 
the strong resolvent sense by a sequence of Hamiltonians with point interac- 
tions supported by a finite set Y C Y . Hence we get a good approximation 
of the spectrum cutting the graph arms to a finite length, large enough. The 
most natural choice of subsets is to consider stars with finite number M 
of point interaction on each arm (and with the central point.) This opera- 
tor if at (a, Y) has the essential spectrum equal to R+ and at most MN + 1 
negative eigenvalues, due to the presence of point interactions. The lower 
eigenvalues, those smaller than E , converge to the eigenvalues of H^{a,[3) 
as M increases, while the rest approximates the negative part of the essential 
spectrum of H N (a,j3). 

The eigenvalues can be obtained as solution to the implicit equation anal- 
ogous to (|4.2j), where 



{« ; z) = ± («>(!) -ln(f)) 
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G (ai, dj] z) = — K (>/^z \di - dj\) 

for i, j = 1, . . . , (MiV + 1). Once we have an eigenvalue Zq it is easy to write 
the appropriate eigenfunction f{x). From | AGHH , Sec. II. 1] we know 

MN+l 
3=1 



where dj are elements of an eigenvector of A(a,Y; zq) corresponding to zero 
eigenvalue. 



5.2 A broken line, N = 2 

Let us start off with a two-arm "star" . Consider 20 point interactions on each 
arm and I = 1. We have proven above that number of eigenvalues of Hn{c(, (3) 
below a fixed energy value increases as (3 goes to zero. Numerical results for 
H 2 (a = 0,Y) with (3 < tt/10 plotted in Fig. |l] agree with this statement; 
they also hint that all eigenvalues are strictly increasing as functions of (3. 
For larger (3, we have a similar situation, see Fig. [2|. We notice that for 
eigenvalues close to the threshold E the approximation by finite-arm star 
with M = 20 becomes insufficient as the picture shows. The eigenvalues 
above the threshold will approximate the continuous spectrum as M — ► oo. 




0.1 0.15 0.2 0.25 0.3 

s 



Figure 1: Several lower eigenvalues of H 2 (a = 0,Y) for small (3. The dotted 
line is the threshold E . 
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Figure 2: The dependence of several lower eigenvalues of H 2 (a = 0,Y) on 
the angle (3. The dotted line is the threshold E . 



P=2/3 it 





10 20 




Figure 3: Eigenfunctions of six lower states of H 2 {a = 0,Y) for (3 = 2/3n. 
Only the ground state has energy below the threshold E . The bold curves 
represent the nodal lines, the contours showing horizontal cuts correspond to 
a logarithmic scale. 
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P=Jt/20 



0^: 



10 15 20 





cj(o 




2 








10 15 20 



10 15 20 




Figure 4: Eigenfunctions of ten lowest states of if2(ct = 0, Y) for /3 = n/20. 
First five states correspond to eigenvalues of H 2 (0,/3), the rest would belong 
to the essential spectrum in the limit M — > oo. The bold curves represent the 
nodal lines, the contours showing horizontal cuts correspond to a logarithmic 
scale. 



In the proof of Theorem [4.1| we observed that pairs of mutually close 
point interactions are crucial for the lower part of discrete spectrum if f3 is 
small. This behaviour can be also demonstrated on the corresponding eigen- 
functions, compare the contour graphs in Fig. ^ to those in Fig. |4]. They 
represent several lowest states of H 2 (a = 0,Y) for the angles (3 = 2/3-7T and 
(3 = 7r/20, respectively. As indicated above, higher eigenfunction will corre- 
spond to the continuous spectrum in the limit M — » oo. This applies to the 
eigenfunctions in Fig. || except the first one which approximates the ground 
state of -^(O, 2/37r). In this case it is the only state which approximates 
an eigenstate of the infinite star. For a much smaller (3 in Fig. f| there are 
five eigenvalues below the threshold (number 0, 1, 2, 3, 4 in the figure), which 
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Figure 5: Ground state energy of H 3 (a = 0,Y). 

correspond to the discrete spectrum of H 2 (0, 7r/20). Notice that the remain- 
ing eigenfunctions resemble a standing-wave pattern along the graph arms as 
one would expect from an approximation from a generalized eigenfunction. 
It may seem that the graph number 5 in Fig. |] gives rise to a bound state 
too, but this only due to an insufficient length M in our approximation. 

5.3 A three- arm star 

Here we consider 10 point interactions on each arm and we put / = 1 again. 
The behavior of eigenvalues is similar to the two-arm case, but the spectrum 
depends of two parameters (3\ and /?2- The minimum binding is achieved in 
the symmetric the graph of ground state energy of H 3 (a = 0, Y) in 

Fig. |5| shows. We see that the eigenvalue does not change much unless one 
of the angles becomes small. The ground state for the symmetric star, f3\ = 
(32 = 2/37T, is illustrated in Fig. || we see the logarithmic singularities at the 
point-interaction sites and the overall exponential decay of the eigenfuction 
along the graph arms. 
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Figure 6: Ground state of H^(a = 0,Y) for f3\ — 02 — 2/37T which approxi- 
mates the ground state of H 3 (0, (2/37T, 2/37r)). 




Figure 7: Eigenfunction of third excited state of H 10 (a = 0,Y) for 
(3 = (0.3,0.7,1.5,1.8,2.5,3,4,4.4,5.2). It approximates an eigenfunction of 
H w (0,(3). The bold curve represents the nodal line, the contours showing 
horizontal cuts correspond to a logarithmic scale. 
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Figure 8: Eigenfunction of third excited state of H w (a = 0,Y) for (3 = 
(7r/5, 2/57T, . . . , 9/57r). It approximates an eigenfunction of Hi (0,/3). The 
bold curve represents the nodal line, the contours showing horizontal cuts 
correspond to a logarithmic scale. 

5.4 Larger N 

In a similar way one can treat star graphs with larger N. In order not 
to overload the paper with the illustrations, we restrict ourselves to a single 
example with N = 10. The nodal line plots on the above pictures call to mind 
the question whether an eigengunction can have a closed nodal line. Such 
states can be found in spectrum of H w (a, Y) for a = 0, as it is illustrated 
in Fig. [7] and Fig [8] for a non-symmetric and symmetric star. One of many 
mathematical questions which can be asked within the present model is about 
the minimum number N for which this is possible. 
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